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ABSTRACT 

We examine the correlation between globular cluster (GC) color and magnitude using HST/ACS 
imaging for a sample of 79 early-type galaxies (—21.7 < Mb < —15.2 mag) with accurate surface- 
brightness fluctuation distances from the ACS Virgo Cluster Survey. Using the KMM mixture modeling 
algorithm, we find a highly significant correlation, 7 2 = rf ^~ z ^ = —0.037 ± 0.004, between color and 
magnitude for the subpopulation of blue GCs in the co-added GC color-magnitude diagram of the 
three brightest Virgo cluster galaxies (M49, M87 and M60). The sense of the correlation is such 
that brighter GCs are redder than their fainter counterparts. For the single GC systems of M87 and 
M60, we find similar correlations; M49 does not appear to show a significant trend. There is no 
correlation between (g — z) and M z for GCs belonging to the red subpopulation. The correlation 

7 g = dt " 9 dg z ^ for the blue subpopulation is much weaker than j z . Using Monte Carlo simulations, we 
attribute this finding to the fact that the blue subpopulation in M g extends to higher luminosities 
than does the red subpopulation, which biases the KMM fit results. The highly significant correlation 
between color and M z , however, is a real effect: this conclusion is supported by biweight fits to the 
same color distributions. We identify two environmental dependencies which influence the derived 
color-magnitude relation: (1) the slope of the color-magnitude relation decreases in significance with 
decreasing galaxy luminosity, although it remains detectable over the full luminosity range of our 
sample; and (2) the slope is stronger for GC populations located at smaller galactocentric distances. 
These characteristics suggest that the observed trend is, at least partially, shaped by external agents. 
We examine several physical mechanisms that might give rise to the observed color-magnitude relation 
including: (1) presence of contaminants like super-clusters, stripped galactic nuclei, or ultra-compact 
dwarfs; (2) accretion of GCs from low-mass galaxies; (3) stochastic effects; (4) the capture of field 
stars by individual GCs; and (5) GC self-enrichment. Although none of these scenarios offers a fully 
satisfactory explanation of the observations, we conclude that self-enrichment and field-star capture, 
or a combination of these processes, offer the most promising means of explaining our observations. 
Subject headings: galaxies: clusters: individual: Virgo - galaxies: dwarf - galaxies: fundamental 
parameters - galaxies: nuclei - globular clusters: general 
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1. INTRODUCTION 

It is by now an accepted paradigm that the formation 
of globular clusters (GCs) is closely linked to the forma- 
tion of their host galaxies (e.g. Searle & Zinn 1978, Ash- 
man & Zepf 1992, Kissler-Patig 1997, Hilker et al. 1999b, 
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Forbes et al. 1997, Cote et al. 1998, Beasley et al. 2002). 
In recent years, observational data for extragalactic GC 
systems have made a huge leap forward in quality, thanks 
mainly to high-resolution imaging from the Hubble Space 
telescope, and to wide field imaging and multi-object 
spectroscopy from 8m-class ground based telescopes (see 
West et al. 2004 and references therein). 

As a rule, studies of extragalactic GCs have tended 
to focus on correlations between the global photomet- 
ric properties of the GC systems (GCSs) and t hose o f 
their host galaxy. For instance, Peng et al. (2006), 
examined the color distributions for GCs belonging to 
100 galaxies observed in the ACS Virgo Cluster Survey 
(ACSVCS; Cote et al. 2004; see also below), finding 
nearly all of these galaxies (— 22 < Mb < —15 mag) to 
possess bimodal, or at least asymmetric, GC color dis- 
tributions. Consistent with previous findings (e.g., Geb- 
hardt & Kissler-Patig 1999; Kundu et al. 2001; Larsen 
et al. 2001), the mean GC color was found to correlate 
with galaxy luminosity in the same sense as the well- 
known color-magnitude relation for early type galaxies 
(e.g. Bower, Lucey & Ellis 1992; Karick et al. 2003; 
Hilker et al. 2003; Ferrarese et al. 2006). A cor- 
responding GC color - galaxy color relation was also 
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found by Peng et al. l)2006t) . confirming previous findings 
(e.g., Burgarella et al. 2001; Larsen et al. 2001; Lotz 
et al. 2004). 

At the same time, analyses of the color magnitude 
diagrams for the GCSs of early-type galaxies revealed 
some tentative evidence for a dependence of ind ividual 
GC color on luminosity. Ostrov et al. l)1998() found 
from Washington photometry of NGC 1399 that the color 
peaks of the GC color distribution seem to merge at the 
highest luminosities. Dirsch et al. l)2003t) confirmed this 
finding for the same galaxy based on an expanded sample 
of GCs. There have also been hints of such a behavior 
in other galaxies such as NGC 5128 (Harris et al. 2004) 
and M87 (Whitmore et al. 1995; Larsen et al. 2001). 

With the improved imaging capabilities of HST pro- 
vided by the Advanced Camera for Surveys (ACS; Ford 
et al. 1998), it has become poss ible to investigate this 
issue in more detail. Harris et al. i|2006f) used ACS imag- 
ing to investigate the distribution of GCs in the (I,B — I) 
color magnitude plane for eight brightest cluster galaxies 
in the range 1800 < cz < 3200 km s _1 . For the four 
galaxies with the largest number of GCs — and also the 
joint sample of eight GC systems — they found a signif- 
icant trend towards redder colors for increasing luminos- 
ity for the subpopulation of blue GCs. This finding led 
them to suggest a mass-metallicity relation of the form 
Z oc M ' 55 for this subpopulation. Curiously, no color- 
magnitude relati on was found for the red GCs. 

Strader et al. (2006) used the publicly available data 
from the ACSVCS to look for color-magnitude trends in 
the three brightest Virgo Cluster galaxies M49, M87 and 
M60. For the blue GC populations of M87 and M60, they 
found a significant trend of redder color with brighter 
GC z-band luminosity similar to that found by Harris 
et al. (200Q). They too interpreted this trend as evidence 
for a mass-metallicity relation among the blue (metal- 
poor) GCs and suggested self-enrichment as a possible 
explanation. No significant color-magnitude trend was 
found for M49, the brightest member of the Virgo Clus- 
ter. 

Both the Harris et al. (|2006t) and Strader et al. ( 2006) 
studies used heteroscedastic mixture models to fit the 
colors of an assumed double-Gaussian color distribution 
to the observed GC color-magnitude diagrams (CMDs) 
as a function of magnitude. Harris et al. use KMM (e.g. 
Ashman, Bird & Zepf 1994), while Strader et al. use Nmix 
(Richardson & Green 1997). Two important issues that 
were not addressed by these studies are the subject of 
this paper: first, we will investigate the color-magnitude 
correlations not only in the red (i.e., z), but also in the 
blue band (i.e., g); second, we will investigate how the 
blue peak slope varies as a function of host galaxy lumi- 
nosity and g alactocentric distance. Note that the Harris 
et al. ( 20QiJ study targeted only brightest cluster galax- 
ies, so it spanned a relatively limited range in galaxy 
luminosity. 

We finally point out that there is recent evidence 
for the existence of a GC color-magnitude relation also 
in sp iral g alaxies. Based on ACS imaging, Spitlcr 
et al. l)2006l) found a significant color-magnitude trend for 
the Sombrero galaxy (NGC 4594) implying Z oc M 30 , 
which is a so mewh at weaker trend tha n sug gested by 
Harris et al. l|2006|) and Strader et al. (2006) for giant 
ellipticals. 



The aim of this paper is to investigate the color- 
magnitude trend of GCs over the entire host galaxy lumi- 
nosity range covered by the ACSVCS. Special emphasis 
will be put on environmental dependencies of this trend 
and checks of the fitting technique. The paper is struc- 
tured as follows: In <J21wc outline how GC candidates are 
selected in the ACSVCS. In Sglwe compare KMM fits to 
simulated CMDs with fits to observed CMDs. In 21 we 
investigate the environmental dependence of the color- 
magnitude trend. In 3Hl we examine a variety of physical 
effects that may create these trends. We finish this paper 
with the summary and conclusions in ^1 

2. SELECTION OF GLOBULAR CLUSTERS FROM THE 
ACSVCS 

The ACSVCS represents the most complete and ho- 
mogeneous study of extragalactic GCs performed to 
date. The sample consists of 100 early-type (E, SO, 
dE, dE,N, dSO) members of the Virgo cluster. Each 
galaxy was imaged in the F475W and F850LP filters 
(~ 5475 an d z§5o, respectively) for a total of 750 sec- 
onds and 1210 seconds, respectively. This filter combina- 
tion gives roughly a factor-of-two improvement in wave- 
length baseline and metallicity sensitivity compared to 
the "canonical" (V — I) color index. The identification 
of bona-fide GCs from these images is performed in th e 
size-magnitude plane as described in Peng et al. (2006). 
This selection procedure — which is possible because the 
half-li ght ra dii of GCs are marginally resolved (Jordan 
et al. i|2005|n at th e distance of Virgo (d = 16.5 Mpc; 
Tonry et al. 1)200 lft : Mei et al. l|2005l^ — greatly re- 
duces contamination from foreground stars and back- 
ground galaxies. 

The ACS images had been reduced using a dedicated 
pipeline described in detail by Jordan et al. (2004ab; Pa- 
pers II and III). In brief, the reduction steps are: im- 
age combination, galaxy modeling and model subtrac- 
tion, rejection of obvious background galaxies, and the 
measurement of magnitudes and sizes for candidate GCs 
using the program KINGPHOT (Jordan et al. 2005; Pa- 
per X). The result is a catalog of integrated g and z 
magnitudes, (g — z) colors and half-light radii, r/j, for 
each candidate GC. Magnitudes and colors are corrected 
for foreground extinction using the reddening maps of 
Schlegel et al. lTl99l . 

To estimate the contamination by background galaxies, 
an identical reduction procedure was applied to 17 blank 
high-latitude fields observed by ACS in the g and z fil- 
ters. Because background galaxies are typically larger at 
comparable luminosities (or, conversely, fainter at com- 
parable sizes) , a statistical decontamination is performed 
in size-magnitude space (see Peng et al. 2006 and Jordan 
et al. 2006, in prep, for details): to this end, first a non- 
parametric density model is fitted to the size-magnitude 
distribution of sources detected in background fields. For 
the GC candidate magnitude distribution, a Gaussian lu- 
minosity function is assumed. The GC candidate size 
distribution is fitted by a non-parametric kernel esti- 
mate plus a power law tail. The resulting surface densi- 
ties in the size-magnitude plane for background sources 
and GCs then make it possible to assign a GC probabil- 
ity, Pgc, to every source detected in the galaxy image. 
All candidates with probabilities V %c > 0.5 are consid- 
ered GCs. The remaining contamination by background 
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Fig. 1— Color-magnitude diagram for GCs in M87 (= NGC4486 
= VCC1316) using M z (left panel) and M g (right panel). Dashed 
lines at (g—z) = 1.0 and (g—z) = 1.4 are shown as a visual aid; they 
correspond to the approximate median colors of the red and blue 
subpopulations. Circles and asterisks indicate KMM fitting results. 
Circles show KMM results for initial guesses of 0.8 and 1.2 mag for 
the two peaks; asterisks show the results obtained for initial guesses 
of 1.0 and 1.4 mag. 

galaxies that are assigned V gc > 0.5 is negligible for the 
sake of this study. We have tested this by statistically 
cleaning the GC CMDs for each galaxy, using the color- 
magnitude distribution of sources with V gc > 0.5 in the 
17 blank background fields. 

3. ANALYSIS OF THE COLOR MAGNITUDE DIAGRAMS 
3.1. Results from KMM 

Figure □ shows M g -(g - z) and M z -(g - z) CMDs for 
the GCS of M87 (= NGC4486 = VCC1316), the cen- 
tral galaxy of the Virgo Cluster and the galaxy with the 
singlest largest population of GCs. For comparison, Fig- 
ure |2 shows the co-added CMDs for the 79 ACSVCS 
galaxies that have reliable SBF distances (Mei et al. 2005, 
2006) and which have a discernible GC population (see 
Peng et al. 2006). We use the SBF distance moduli from 
Mei et al. (2005, 2006), which have an internal precision 
of w 0.07 mag, to assign absolute magnitudes to the asso- 
ciated GC candidates. To better visualize trends across 
host galaxy magnitude, the 79 galaxies are subdivided 
into four broad luminosity groups such that the number 
of GCs per group falls in the range N gc ~ 2400 to 3300. 

Luminosity Group 1 consists of the three bright- 
est AVCVCS galaxies M49 (= NGC4472 = VCC1226; 
M B = -21.7 mag), M87 (= NGC4486 = VCC1316; 
M B = -21.4 mag), and M60 (= NGC4649 = VCC1978; 
M B = -21.2 mag). Group 2 consists of six ACSVCS 
galaxies with —21 < Mb < —20 mag. Group 3 consists 
of 17 ACSVCS galaxies with -20 < M B < -18.4 mag. 
The remaining 53 ACSVCS galaxies with -18.4 < M B < 
— 15.2 mag make up Group 4. The three galaxies in 
Group 1 are all so-called "core-Sersic galaxies" (Ferrarese 
et al. 2006), while those in Groups 2 and 3 are a mixture 
of core-Sersic and Sersic elliptical and lenticular galaxies. 
The 53 galaxies in Group 4 consist mainly, but not ex- 
clusively, of early-type dwarfs (e.g., dE, dE,N, dSO and 
dS0,N). 



We have applied KMM to the combined CMDs for these 
four groups, as well as to the individual CMDs for the 
three galaxies in Group 1. To this end, we subdivided the 
CMDs into luminosity bins, each bin containing the same 
number of data points. For the CMDs of M49 and M60, 
the size of each luminosity bin is 100, for M87 the size is 
150, and for the combined CMDs the bin size is 200 GCs. 
KMM was then run on the GC color distribution i n eac h 
luminosity bin. As in the paper by Harris et al. (2006), 
we use the heteroscedastic option of KMM: we allow for 
Gaussians of different widths to be fitted to the blue 
and red GC populations. The fitted mean positions are 
plotted over the respective CMDs in Figures ^ an d EJ 
We show the KMM results for two different initial guesses 
for the mean colors of the blue and red GC populations: 
one for the pair (g — z) = 0.8 and 1.2 mag, another one 
for (g — z) = 1.0 and 1.4 mag. These different guesses 
were chosen to cover: (1) the range of peak positions 
between CMDs of the various luminosity groups; and 
(2) the range caused by the color-magnitude trend itself 
(see Figure . Linear slopes between color and absolute 
magnitude are obtained from ordinary least square fits to 
the KMM peak positions for both pairs of initial guesses, 
fitting color as a function of magnitude. GCs fainter than 
M z = — 7.7 mag and M g = — 6.7 mag were excluded from 
the fits, to avoid effects of photometric incompleteness on 
the measured slopes. The mean of both fits is adopted 
as the final best-fit slope. The fit errors of the slope are 
derived from re-sampling the points using as dispersion 
the observed scatter around the fitted relation. We note 
that this scatter is for all luminosity groups larger (by 
30% for luminosity group 1 and by 100% for group 4) 
than the average statistical error estimate for the peak 
position, given by Here, a is the KMM width of the 

Gaussian peak and N the average number of GCs per 
magnitude bin associated to the peak. We are therefore 
confident that our error estimates are not too optimistic. 
For the few cases where the difference between the results 
obtained using the two initial guesses was larger than the 
formal fit errors, we adopted the difference as the error. 
The resulting slopes of those fits are shown in Table 
We first consider M z vs. (g — z) for Group 1. In this 

case, we detect a very significant slope "/ z = = 
—0.037 ± 0.004 for the blue peak position. When re- 
stricting the fit to M z > —10 mag (i.e., rejecting the 
two brightest magnitude bins in Figure EJ, the result is 
j z = -0.033 ± 0.007. This is significant at the 4.5cr 
level, indicating that the color-magnitude trend is not 
restricted to just the brightest GCs. For M87 and M60 
alone, a comparable slope is derived, albeit at some- 
what lower (w 3a) significance. The slope for M49 is 
consistent with zero , in agreement with the findings of 
Strader et al. ( 2006). When excluding M49 from Group 
1, the slope for the blue population rises slightly to 
j z —0.040. The non-detection of a color-magnitude 
trend for M49 is intriguing, since this galaxy has a simi- 
lar number of GCs as M60 and is similar to M60 and M87 
in terms of luminosity, color, surface brightness and GC 
dynamics (Peng et al. 2006; Ferrarese et al. 2006; Mei 
et al. 2005; Cote et al. 2001; 2003). The fact that M49 
is located at the center of its own sub-cluster — being 
offset by « 4.5° from the Virgo Cluster center — may 
hint at some environmental effect on the strength of the 
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FlG. 2. — The same as Figure [Tl except for the combined sample of GCs in each of the four host-galaxy luminosity groups (see Tabic E). 
(Upper left) The co-added CMD for GCs belonging to the three brightest ACSVCS galaxies, having Mb < —21 mag (Group 1). (Upper 
right) The co-added CMD for GCs belonging to the six ACSVCS galaxies with -21 < M B < -20 mag (Group 2). (Lower left) The 
co-added CMD for GCs belonging to the 16 ACSVCS galaxies with -20 < M B < -18.4 mag (Group 3). (Lower right) The co-added CMD 
for GCs belonging to 53 ACSVCS galaxies with —18.4 < Mb < —15.2 mag (Group 4). The luminosity groups have been chosen such that 
the number of GCs in each co-added CMD are comparable (i.e., N gc = 2400 — 3300). 



observed color-magnitude trend. 

In contrast to the strong correlation ^ z exhibited by the 
blue GCs in Group 1, j z is insignificant for the red GCs 
in this group. However, the situation changes when using 
M g instead of M z as magnitude in the color-magnitude 

diagram (see Table The slope 7 9 = d< " 9 dg z " > for the 
blue population is more than 50% smaller than "f z and, 
as a result, has a lower level of significance (1.7a.) Con- 
sidered individually, j g of the blue population becomes 
insignificant for all three galaxies in this group. In turn, 
the slope j g of the red population becomes marginally 
positive — i.e., opposite in sense to that defined by j z 
for the blue population. 

3.2. Monte Carlo Simulations to test KMM 

The difference in slopes between z and g raises the 
question as to what the "true" slope actually is. To in- 



vestigate this issue in more detail, we compare in FigureEl 
the actual measurements for Group 1 with two sets of 20 
Monte-Carlo simulations for the M87 CMD. The first set 
of CMDs has an implemented slope 7 2 = —0.040 for the 
blue GCs 10 and 7 Z = for the red ones; the other set has 
7 2 = implemented for cither GC population. These ar- 
tificial CMDs were created by assuming the average peak 
widths derived by KMM for the blue and red GCs, taking 
into account also the photometric errors on the colors. 
The luminosity function is sampled from a kernel esti- 
mator of the original CMD of M87, separately for each 
color peak. Figure 01 compares the real and simulated 
M87 CMD with an implemented slope. Note an impor- 
tant feature of those CMDs: the red peak is populated to 

10 Note that the slope 7 9 corresponding to 7 Z = —0.040 is of 
course not identical to 7 Z , but this difference is negligible for the 
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Fig. 3. — KMM result for the co-added CMD of the three brightest ACSVCS galaxies (Group 1, see Figure[3]and Table0 compared to 
KMM results for simulated M87 CMDs. Filled circles indicate KMM fitting results for the real data with initial guesses of 1.0 and 1.4 mag for 
peak positions (see Figures [Tl and |21 . Dots indicate KMM peak fittings for 20 simulated CMDs. (Left Panel) Simulated CMDs with a slope 

of ■jz = d ^ 9 dz z ^ = —0.04 implemented for the blue subpopulation; no correlation has been implemented for the red subpopulation. Note 
the artificial bending of the red sequence toward blue colors at bright luminosities, especially for M g . (Right Panel) Simulated CMD with 
no color-magnitude relation implemented for either population. 
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Fig. 4. — Comparison between real (left panel) and simulated (right panel) CMD for GCs in M87. The simulated CMD has a slope of 
7 Z = —0.04 input for the blue subpopulation. 



brighter absolute luminosities in M z than in M g , while 
the blue peak does not show such a difference. This is 
because, relative to GCs in the blue population, GCs in 
the red population are brighter in M z than in M g . 

Table Ogives the slopes fitted to the mean colors from 
KMM for the two simulated sets of M87 CMDs, together 
with the slopes for the real data. This comparison shows 
that both the CMD of M87 alone, and that for the co- 
added CMD of Group 1, are inconsistent with a magni- 
tude independent color for the blue GCs. The data are 
consistent with a simulated slope of -0.040 for the blue 
GCs and no slope for the red GCs. 

Note the positive slope found for the red GCs in both 
the real and simulated data in M g . Averaging over the 



20 simulated CMDs, the mean red slope of 7 g ~ +0.03 is 
different from zero at the 5.7a level. That finding is note- 
worthy since the simulations do not implement a slope 
for the red peak. This exercise therefore demonstrates 
an important caveat of the KMM fitting: at a given mag- 
nitude, as one of the two peaks becomes de-populated 
(in this case the red peak in M g ), its KMM position will 
artificially be dragged towards the position of the more 
prominent peak. As a consequence, the mean color of the 
blue population will be artificially shifted towards bluer 
colors. This results in a weaker slope for the blue peak 
in M g than in M z . In M z , both peaks are populated up 
to comparable absolute luminosities, and therefore the 
fitted slopes for both red and blue peaks correspond to 
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TABLE 1 

Color-magnitude trends for red and blue globular clusters in ACS Virgo Cluster Survey galaxies, determined with 

kmm FITS 



Sample 

M49 (M B = -21.7) 
M87 (M B = -21.4) 
M60 (Mb = -21.2) 



-0.008 
-0.042 
-0.028 



7z,bluc 

± 0.024 
± 0.015 
± 0.009 



0.013 
0.003 
0.012 



,rcd 

± 0.018 
± 0.023 
± 0.018 



0.021 
•0.016 
0.003 



7g,bluc 

± 0.018 
± 0.023 
± 0.033 



0.043 
0.025 
0.072 



7g,red 

± 0.008 
± 0.017 
± 0.017 



Simulated M87 NoCM 
Simulated M87 CM 



0.002 
-0.039 



± 0.012 
± 0.009 



0.010 
0.010 



± 0.015 
± 0.014 



0.008 
0.028 



± 0.015 
± 0.010 



0.028 
0.031 



± 0.022 
± 0.024 



Groupl (-21.7 < M B < -21) 
Group2 (-21 < M B < -20) 
Group3 (-20 < Mb < -18.4) 
Group4 (-18.4 < Mb < -15.2) 
Groupl, inner regions 
Groupl, outer regions 
Group 1+2, inner regions 
Groupl+2, outer regions 



-0.037 
-0.033 
-0.032 
-0.009 
-0.061 
-0.026 
-0.066 
-0.023 



± 0.004 
± 0.011 
± 0.012 
± 0.013 
± 0.014 
± 0.008 
± 0.007 
± 0.013 



0.001 
0.001 
-0.043 
0.028 
-0.006 
0.003 
-0.024 
0.009 



± 0.007 
± 0.014 
± 0.024 
± 0.028 
± 0.014 
± 0.020 
± 0.006 
± 0.013 



0.015 
■0.022 
0.010 
0.025 
-0.007 
0.004 
-0.022 
0.018 



± 0.009 
± 0.016 
± 0.007 
± 0.012 
± 0.018 
± 0.015 
± 0.020 
± 0.007 



0.032 
0.042 
0.055 
0.058 
0.059 
0.050 
0.043 
0.029 



± 0.019 
± 0.018 
± 0.029 
± 0.040 
± 0.018 
± 0.014 
± 0.018 
± 0.013 



Notes: For the samples in column 1, columns 2 to 5 give the slopes 7 between (g — z) and the respective magnitude of the blue 
and red GC subpopulations, as derived from KMM determined peak positions. Errors come from random resampling of the data 
points using their measured dispersion around the fit. For the two simulated sets of CMDs, "CM" refers to an implemented 
slope 7 Z = —0.040 in the blue peak and no slope in the red peak. "NoCM" refers to no implemented slope. Errors quoted 
for the simulations are rms around the mean for 20 simulated M87 CMDs. The faint limiting magnitude for fitting was —7.7 
mag in M z and —6.7 mag in M g . The separation between inner and outer region for the lower four samples was defined at a 
galactocentric distance of 65". By that, the inner and outer sample were of approximately equal size. 



TABLE 2 

Color-magnitude trends for red and blue globular clusters in ACS Virgo Cluster Survey galaxies, determined with 

biweight fits 



Sample 

M49 (M B = -21.7) 
M87 (M B = -21.4) 
M60 (M B = -21.2) 



-0.004 
-0.031 
-0.024 



7z,bluc 

± 0.002 
± 0.007 
± 0.006 



0.003 
0.013 
0.001 



7z,rcd 

± 0.003 
± 0.005 
± 0.005 



0.005 
-0.031 
0.027 



7g.bluc 

± 0.006 
± 0.013 
± 0.007 



0.010 
0.017 
0.017 



± 0.004 
± 0.006 
± 0.002 



Simulated M87 NoCM 
Simulated M87 CM 



0.008 
-0.029 



± 0.007 
± 0.006 



0.001 
0.013 



± 0.005 
± 0.004 



0.008 
0.026 



± 0.004 
± 0.007 



0.009 
0.022 



± 0.007 
± 0.003 



Groupl (-21.7 < Mb < -21) 
Group2 (-21 < M B < -20) 
Group3 (-20 < Mb < -18.4) 
Group4 (-18.4 < M B < -15.2) 
Groupl, inner regions 
Groupl, outer regions 
Groupl+2, inner regions 
Groupl+2, outer regions 



-0.025 
-0.025 
-0.020 
-0.019 
-0.035 
-0.022 
-0.032 
-0.025 



± 0.002 
± 0.004 
± 0.002 
± 0.004 
± 0.005 
± 0.004 
± 0.004 
± 0.005 



0.009 
0.017 
-0.003 
0.011 
0.012 
0.003 
0.016 
0.010 



± 0.003 
± 0.005 
± 0.009 
± 0.007 
± 0.004 
± 0.009 
± 0.002 
± 0.006 



0.027 
0.023 
-0.001 
-0.021 
-0.043 
-0.020 
-0.029 
0.022 



± 0.004 
± 0.003 
± 0.005 
± 0.004 
± 0.006 
± 0.005 
± 0.006 
± 0.003 



0.021 
0.025 
0.018 
0.012 
0.022 
0.016 
0.026 
0.017 



± 0.004 
± 0.006 
± 0.010 
± 0.012 
± 0.010 
± 0.004 
± 0.004 
± 0.007 



Notes: Analogous to Table Q Difference is that peak positions are now defined by the biweight of the color distribution blue- 
and redwards of a luminosity independent limiting color, see Fig. |K| Fig. |S| graphically compares the fits to KMM and biweight 
peak positions for all four luminosity groups. 



the simulated ones. That is to say, only the weaker slope 
in M g is caused by a KMM bias. KMM does not artificially 
introduce the very significant negative slope in the blue 
M z peak: this is clearly a real trend that is intrinsic to 
the data. 

When interpreting the slope as a mass-metallicity re- 
lation for GCs in the blue peak, we derive a scaling 
Z cx M°- 48±0 08 by combining the slope j z =-0.037 found 



in the previous section with 
d[Fe/H] 



5.14 dex mag 



from Peng et al. |2QQf|) and assuming a constant mass- 
to-light ratio over the luminosity range investigated. 
This is in good agreement with the scaling relation 
Z oc A/ 55 found by Harris et al. l|200^ . 

As an alternative to KMM, we have calculated the bi- 
weight location (Beers et al. 1990) for the GC subpopula- 
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tions blue-wards and red-wards of a magnitude indepen- 
dent limiting color. This limiting color is the mean of the 
KMM position of the blue and red populations, weighted 
by the average width of either peak. The results for this 
alternative method are shown in TableEl As an example, 
we show in Fig. the peak positions for Groups 1 and 4 
in M z . Generally speaking, the errors for this biweight 
method decrease significantly compared to the KMM data 
points, and the results confirm the significant slope found 
for the blue GC population. They also confirm that no 
significant correlation is found for M49. Tableland Fig- 
ure|S|show that the difference between fits in M g and M z 
is considerably smaller than in the case of KMM. 

Nevertheless, this alternative method is itself not im- 
mune to biases, as is clear when it is applied it to the 
simulated CMDs (see Table Eland Figure EJ). In the case 
of the blue population exhibiting a negative slope, a fixed 
limiting magnitude between both populations will sys- 
tematically decrease the absolute slope of the blue pop- 
ulation, and introduce an artificial positive slope for the 
red. 

4. ENVIRONMENTAL EFFECTS 

Having examined how specific features of the CMD 
may influence the fitted color-magnitude relations for the 
blue and red GC populations, we now investigate the ex- 
tent to which possible environmental effects may influ- 
ence the observed trends. 

4.1. Dependence on Host Galaxy Luminosity 

Figure shows CMDs for Groups 1-4 with the 
KMM peaks over-plotted. Derived slopes for the color- 
magnitude relations in each group are presented in Ta- 
bles and [21 The dependence of slope on host galaxy 
magnitude is illustrated in Figure El 

The KMM fits in M z show that the slope of the blue 
GCs in Groups 2 and 3 are indistinguishable from that 
for Group 1; in all three groups, the slope is inconsistent 
with zero at about the 3cr level. This indicates that the 
color-magnitude trend among metal-poor GCs is not re- 
stricted to just the brightest galaxies. Only for Group 4 
does the measured slope of the blue population become 
consistent with zero. The biweight fits in the lower pan- 
els of FigureElpaint a similar picture: the slope decreases 
only gradually from Group 1 to 4, and apparently does 
not drop to zero. Given that the biweight fits systemat- 
ically underestimate the blue slope, this is an indication 
that the color magnitude trend persists over the entire 
magnitude range of the ACSVCS galaxies. 

The trend for the red GCs determined using KMM be- 
comes progressively more uncertain for fainter galaxy 
magnitudes. This is because of the decreasing fraction of 
red (metal-rich) GCs with decreasing host galaxy lumi- 
nosity (Peng et al. 2006). The firmest result from Fig- 
ure El regarding the color-magnitude relation for red GCs 
is that there is no strong M z slope for the two bright- 
est groups; the slopes for the two faintest galaxy groups 
have large uncertainties. The large positive slopes found 
in M g are particularly prone to KMM biases, a conclusion 
that is confirmed by looking at the more robust biweight 
fit results. The comparison between the biweight fits to 
the real data and the simulations shows that the slope 
of the red population is consistent with zero over the full 
luminosity range. 



4.2. Dependence on Galactocentric Distance 

An interesting feature of the observed color-magnitude 
relations is illustrated in Figured Here we subdivide the 
GC samples of Groups 1 and 2 by galactocentric radius 
into equally-sized, "inner" and "outer" subsamples. As 
the dividing radius, we choose R = 65", which translates 
to a physical distance of ss 5.2 kpc at the distance of the 
Virgo Cluster. In terms of the effective radii, R e , of the 
host galaxies in Groups 1 and 2, this fixed dividing ra- 
dius corresponds to « 0.3-0.6 R e (Ferrarese et al. 2006). 
Tables Q and report the slopes for the different sub- 
samples. For Group 1, the slope found for the inner GCs 
is a factor of ps 2-3 higher than that for the outer GCs. 
Both the KMM and biweight methods support this conclu- 
sion, at significance levels of 2.2cr and 2.5cr, respectively. 
For the combined Group 1+2 samples, the KMM fits still 
obtain a difference of about a factor of three (at 2.9a sig- 
nificance), while the differences are almost insignificant 
for the biweight fits. 

In Figure 00 we plot color vs. half-light radius for GCs 
with —11.5 < M z < —10 mag, separated into inner and 
outer samples. There is a small size offset between the 
inner and outer sample, in the sense that outer GCs are 
larger by about 10%. This dependence of r^ on galac- 
tocentric radius has already been investigated by Jordan 
et al. (2005). In Figure El we plot the color histograms 
for two luminosity ranges: in the left panel, we show GCs 
in the same luminosity range as in Figure |H1 in the right 
panel, we show GCs in the range — 10 < M z < —9 mag. 
Clearly, the gap between the two peaks is larger among 
the fainter sample of GCs, a consequence of the color- 
magnitude trend. Furthermore, for the brighter GCs, the 
inner and outer color histograms appear to differ more 
than for the fainter GCs. The KS probability with which 
the inner and outer sample of the bright GCs have the 
same underlying distribution is 1%, while it is 14% for 
the faint sample. Indeed, while the bimodality feature 
is still notable for the outer bright GCs, it is slightly 
smeared out for the bright, inner GCs. This suggests 
that the dependence of the color-magnitude relation on 
galactocentric distance is driven mainly by the bright 
GCs (M z < -10 mag). 

FigureElalso shows that the fraction of blue-to- red GCs 
decreases slightly for smaller galactocentric distance (see 
also, e.g., Dirsch et al. 2003). The corresponding frac- 
tions of blue-to-red GCs are 0.45:0.55 and 0.50:0.50 for 
the inner and outer GCs, respectively. It is apparent from 
the color histograms in Figure that this small effect 
cannot be responsible for the steeper color-magnitude 
relation found for the inner sample of GCs. 

4.3. Comparison to CMDs of Other Stellar Systems 

Before proceeding, we pause to compare the distribu- 
tion of GCs in the CMD with those of other stellar sys- 
tems: i.e., the nuclei of ACSVCS galaxies, the ACSVCS 
galaxies themselves, and some candidate ultra-compact 
dwarf galaxies (UCDs) identified in the course of the 
ACSVCS. 

Figure ^3 compares the combined CMD of GCs 
(red circles), compact stellar nuclei (asterisks; Cote 
et al. 2006), ACSVCS galaxies (green circles; Fer- 
rarese et al. 2006) and UCDs (open squares; Hasegan 
et al. 2005). In this figure, we have converted from M z 
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Fig. 5. — (Left Panel) Color-magnitude diagram in M z for GCs in Group 1. Filled circles indicate the biweight estimates for the mean 
GC color bluewards and redwards of the dividing color (g — z) = 1.155, which is indicated as a dotted vertical tick mark. The solid lines are 
least-squares fits to the circles. For comparison, the dashed lines give the linear fit obtained from the KMM colors (Table HI . (Right Panel) 
Same as in the left panel, but for Group 4. The offset between the KMM and biweight fits for the red subpopulation is due to its asymmetric 
color distribution (i.e., an extended tail towards red colors). 

tion, nor with that of the nuclei. Moreover, the color- 
magnitude relations of the nuclei and the galaxies are 
steeper than that of the GCs. Only the faintest nuclei 
appear to be roughly consistent with the color-magnitude 
relation observed for the blue GCs. In any case, we point 
out that the existence of a color-magnitude trend for the 
blue GCs results in a paucity of GCs at high luminosities 
and blue colors: i.e., M z < — 10 mag and (g — z) < 0.9 
mag. It is difficult to understand how contamination of 
the observed CMDs by some unrecognized stellar system 
— stripped nuclei, ultra-compact dwarfs (UCDs), or oth- 
erwise — can lead to an underpopulation of this region 
of the CMD. 
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Fig. 6. — Slopes 7 of the GC color-magnitude relation for the four 
host-galaxy luminosity groups plotted against the respective mean 
host galaxy magnitude. Filled circles and asterisks show results 
for 7 Z and y g , respectively. (Upper Panels) KMM fits (see Table ITt . 
(Lower Panels): Biweight fits (see Table Izl . (Left Panels) Slopes 
for the blue subpopulations. (Right Panels) Slopes for the red 
subpopulations. The horizontal lines indicate the respective fit 
results to the artificial CMD of M87, with a slope of ^ z = -0.04 
implemented for the blue subpopulations and zero slope for the red 
subpopulations. Dashed and dotted lines correspond to -y z and 7 9 , 
respectively. 



to stellar masses usi ng the population synthesis models of 
Bruzual & Chariot l|2003l) under the assumption of a 10- 
Gyr old population and a Chabrier IMF (Chabrier 2003). 
For comparison, we also show the color-magnitude rela- 
tion for early- type g alaxies in the Fornax Cluster from 
Hilker et al. (2003J), after converting from (V — J) t o 
(g—z) using the transformation of Fukugita et al. ( 1996). 

The majority of the blue GCs are not consistent with 
an extrapolation of the galaxy color-magnitude rela- 



5. DISCUSSION 

As a starting point for this discussion, we summarize 
the two main findings thus far, both of which concern 
dependencies of the measured color-magnitude slope on 
external factors: (1) The slope for the blue GCs de- 
creases gradually towards fainter host galaxy luminosi- 
ties, whereas there is no significant color-magnitude re- 
lation for the red GCs; and (2) the slope is steeper for 
GCs which lie closer to the centers of their host galaxies. 

We now discuss a number of scenarios that that might 
help explain the observed trend: (1) contamination by 
super star clusters, stripped nuclei, or UCDs; (2) accre- 
tion of GCs from low-mass galaxies; (3) stochastic effects; 
(4) capture of field stars by individual GCs; and (5) GC 
self enrichment. In fl5.fi! we briefly discuss how the ob- 
served color-magnitude trends appear after transforming 
into the metallicity-mass plane. 

5.1. Super Star Clusters, Stripped Nuclei, and UCDs 

One idea that may explain the slope difference be- 
tween inner and outer GCs is the presence of overlu- 
minous star clusters of intermediate color that may lie 
close to the galaxy centers. Such a putative popula- 
tion would populate the region between the two color 
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Fig. 7. — Radial dependence of the color-magnitude for the combined GC systems of Groups 1 and 2. (Left Panel) GCs having 
galactocentric distances R < 65" (corresponding to d < 5.2 kpc). (Right Panel) GCs having galactocentric distance R > 65". Note the 
stronger trend for the blue subpopulation inside 65" (see also Table HI . 
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Fig. 8. — Half-light radius, r%, plotted against (g — z) for GCs 
belonging to Group 1 and 2, and with magnitudes in the range 
— 11.5 < M z < —10 mag (see also Figure 171. Open circles are GCs 
with galactocentric distances R > 65", filled circles are GCs with 
R < 65". The mean size of the outer GCs is larger than that of 
the inner GCs by about 10% (see also Jordan et al. 2005). The 
cumulative r h distributions of inner and outer GCs are inconsistent 
with each other at the 99.99% confidence level according to a KS 
test. 

sequences at high luminosities. Merged super star clus- 
ters (e.g., Fellhauer & Kroupa 2002, 2005; Kissler-Patig 
et al. 2006) could, in principle, be those objects. They 
are expected to be more extended than their monolithic 
counterparts of equal mass (e.g., Bekki et al. 2004; Fell- 
hauer & Kroupa 2002; Bastian et al. 2005). However, in 
the size vs. color plot shown in Figure |H1 it is the outer 
GCs that are slightly larger. More importantly, there is 
no offset in the size distribution for the inner GCs with 
intermediate colors (i.e., « (g — z) — 1.20 mag). This 
makes the presence of a large number of merged star 
clusters unlikely. 
Another possibility is the presence of "naked" or tidally 



stripped dE,N nuclei, as suggested by Harris et al. (|2£ 
(see also Bassino et al. 1994). The dE,N nuclei of Virgo 
Cluster galaxies do exhibit a color-magnitude trend and 
have colors comparable to those of the blue GCs (see, 
e.g., Lotz et al. 2 004; Cote et al. 2006; Figure UTty . 
Cote et al. (2006) find a median half-light radius of 
r/j = 3.2±0.6 pc for nuclei with —11.5 < M z < —10 mag. 
This is comparable to the sizes of GCs in this same lumi- 
nosity range (Figure |5J). Stripped nuclei might therefore 
be confused with "normal" blue GCs at the highest lu- 
minosities. 

However, as pointed out in H4.3I the overall distribu- 
tion of GCs with the CMD shows a lack of bright, blue 
objects. Very generally, this dearth is not explained by 
the addition of an "extra" population of contaminators, 
but rather requires an underlying trend for most of the 
objects. Stripped nuclei would have to dominate the pop- 
ulation of blue GCs. This seems unlikely to us for several 
reasons. First, the nuclei luminosity function peaks ~ 2- 
4 mag brighter than that of GCs (see Lotz et al. 2004;, 
Cote et al. 2006, see also Figure ITT)|) . More seriously, 
though, is the shear number of nuclei that would be re- 
quired in this scenario. Observationally, there is an up- 
per limit on the number of naked nuclei brighter than 
M v w —11 mag from Jones et al. (2006). These au- 
thors found nine compact objects - UCDs in their nota- 
tion - with Mb < —10.7 mag in a survey that is almost 
complete within the central degree of the Virgo Cluster 
(excluding the 1-2' surrounding M87). Their nine com- 
pact objects are distributed over an area that is ~ 100 x 
the ACS field of view, in qu alitative agreement with the 
predictions by Bekki et al. ( 20 03aj) for the spatial dis- 
tribution of tidally stripped dE,N nuclei. Therefore, a 
presence of dozens, or possibly hundreds, of naked nuclei 
in our ACS fields for the most luminous galaxies seems 
unlikely. 

A final note regarding the possibility of UCDs as con- 
taminators: the two most discussed formation channels 
for these compact stellar systems are stripped dE,N nu- 
clei and merged stellar super-clusters (e.g. Minniti et 
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Fig. 9. — (Left Panel) GC color distributions for the two samples shown in Figure[S] The dotted histogram is the sample of inner GCs 
with —11.5 < M z < —10 mag; the solid histogram is the sample of outer GCs with —11.5 < M z < —10 mag. (Flight Panel) Same as in 
left panel, except for GCs with —10 < M z < —9 mag. The KS probability for the inner and outer GC samples having the same underlying 
distribution is 1% in the left panel and 14% in the right one. To facilitate the comparison, the histograms for the inner GCs have been 
rescaled by a factor of 1.2x in both plots, so that inner and outer sample have equal sizes. 
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Fig. 10. — Color-magnitude diagram for GCs, compact stellar 
nuclei, early-type galaxies and UCDs in the Virgo Cluster. Small 
filled circles (online red) are GCs belonging to Group 1. The short 
dashed line is the KMM fit to the blue GC subpopulation, extrapo- 
lated to brighter luminos ities. Asterisks are the nuclei of ACSVCS 
galaxies from Cote et al. J2006). Open squares (online blue) are the 
three UCD candidates from Has_egan et al. 1:2005) with ACS imag- 
ing. Large filled circles (online green) are the ACSVCS galaxies 
themselves (Ferrarese et al. 2006). The long dashed line indicates 
the (V — I) color-magnitude relati on for early- type galaxies in the 
Fornax Cluster from Hilker et al. 12003). The Fornax relation was 
transformed into (g — z) accord ing to (g — z) = 2.64(V — I) — 1.72, 
following Fukugita et al. 119961) . The ordinate on the right denotes 
stellar masses, converted from M z using the mean M„/L z ratio at 
metallicities of [Fc/H] = —2.25 and +0.56 dex for a 10- Gyr p opula- 
tion (Bruzual & Chariot 2003), assuming a Chabrier 12003) IMF. 
Due to the low sensitivity of M z on metallicity, M, /L z varies by 
just 40% between the two extreme metallicities, justifying a color 
independent luminosity-mass transformation for illustrative pur- 
poses. 
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Fig. 11. — One test of the accretion scenario, showing a CMD in 
M z and (g — z). Asterisks indicate the KMM fits to the blue and red 
subpopulations for GCs in Group 1 (see Figure Dots show the 
mean GC color as expected from the mean host galaxy luminosity 
in each luminosity bin for the co-added CMDs of Groups 2-4. The 
tilt seen for the dots is the result of two effects: (1) the larger width 
of the GC luminosity function at brighter host galaxy magnitudes; 
and (2) the correlation between the mean GC color and host galaxy 
luminosity. See text for details. 

al. 1998, Hilker et al. 1999a, Drinkwater et al. 2003, 
Ha§egan et al. 2005, Kissler-Patig et al. 2006, Mieske 
et al. 2006). Since these two formation channels cor- 
respond to the two kinds of contaminators discussed in 
this Section, we conclude that also UCDs are unlikely to 
account for the observed trends. 

5.2. Accretion from Low-Mass Galaxies 
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A population of GCs accreted from galaxies with a 
range of masses (e.g. Cote et al. 1998 or Hilkcr 
et al. 1999b) will exhibit a color-magnitude relation if: 
(1) the mean GC color of a single galaxy scales with 
galaxy mass; and (2) more massive GCs preferentially 
form in more massive galaxies. Indeed, both conditions 
are satisfied in actual early-type galaxies (see, e.g., Peng 
et al. 2005; Kundu & Whitmore 2001; Jordan et al. 2006 
in prep, and references therein). Moreover, this simple 
scenario may offer, at least in principal, a way to qual- 
itatively understand the slope differences for individual 
galaxies (e.g., the non-existence of a correlation for M49) 
as a consequence of different accretion/merger histories. 

In order to quantify the tilt in the CMD created by 
the accretion of GCs from low-mass galaxies that do not 
themselves have color-magnitude relations, we consider 
the combined GC sample of Groups 2-4 (see Figure ITT|) . 
We first calculate the mean host galaxy magnitude for 
each bin, and then calculate the color in each GC bin 
according to the relation between host galaxy magnitude 
and mean color of the blue and red GC subpopulations 
(Table 3 of Peng et al. 2006). The outcome of this 
exercise is shown in Figure 1111 As expected, there is a 
(weak) trend between mean color with luminosity for the 
GC subpopulations. 

This trend is compared to the KMM fits for Group 1. 
There is a substantial offset in the sense that the fainter 
luminosity groups have bluer GC colors than Group 1. 
In addition, the simulated slope is substantially smaller 
than measured for the ensemble of galaxies in Group 1. 
All in all, we conclude that the observed color-magnitude 
relations are unlikely to have been created by accretion 
of GCs from low-mass galaxies, although we cannot rule 
out the possibility that accretion has contributed to the 
observed trend. 

5.3. Stochastic effects 

The light of old GCs is in the red optical to near- 
infrared pass-bands dominated by red giant stars. Gen- 
erally, the number of red giants is only a small fraction of 
the total number of stars in a given stellar population. It 
is therefore worth discussing to which extent stochastic 
effects due to small numbers of red giants can skew the 
probability distribution of the integrated GC luminos- 
ity. This may cause a trend of GC color with integrated 
luminosity, because any stochastic effect is weaker for 
higher GC luminosity (larger total number of stars) and 
for bluer bands (lower average stellar luminosity). 

In the works by Cerviho et al. (2002), and Cerviho & 
Luridiana (2004; 2006), the effect of discrete isochrone 
sampling on integrated photometric properties is dis- 
cussed in detail. In Cerviho & Luridiana (2006), the 
width, skewness and kurtosis for the probability distri- 
bution of integrated luminosity of a stellar population 
are calculated as a function of star number and age. The 
faint magnitude limit of our GCs (M z ~ — 7.7 mag) cor- 
responds to about 10 5 M Q , and hence about 2 x 10 5 stars. 
In the /-band - taken to represent the z band due to 
lack of estimates for the latter band -, the photomet- 
ric uncertainty arising from sampling effects amounts to 
about 0.06 mag for such a population (Cerviho & Lurid- 
iana (2006); and Cerviho, private communications). The 
skewness of the distribution in integrated luminosity is 
about 0.25. The product of these two numbers - about 



0.015 mag — is a measure of the difference between mode 
and mean of the probability distribution in integrated 
luminosity. Given that the colour-magnitude trend in 
our data requires an almost ten times larger color shift, 
stochastic effects do probably not have a significant con- 
tribution to the observed trend. 

There is also a phenomenological counterargument 
against the importance of stochastic effects: the color- 
magnitude trend is shaped more by brighter GCs (see 
FigurelHJl, while stochastic effects are stronger for fainter 
GCs. Furthermore, the environmental dependencies of 
the trend cannot be explained by purely intrinsic stochas- 
tic effects. 

5.4. Capture of Field Stars 

Given that the strength of the color-magnitude trend 
depends on local environment, it is reasonable to examine 
the extent to which capture of field stars may help to 
explain the observed slopes. The capture of fiel d stars 
has previously been considered by Bica et al. (|1997f ) 
for the case of a 10 5 M Q star cluster orbiting within the 
Galactic bulge at a distance of 1 kpc. They found that 
such a cluster could absorb a substantial fraction of its 
total mass during its lifetime. Similarly, Kroupa l|1998f) 
has shown that stellar super-clusters created in disks like 
that of the Milky Way can capture several tens of percent 
of their total mass in form of disk stars. 

Field stars of the brightest early-type galaxies are sig- 
nificantly redder than their blue GC population, typi- 
cally by about 0.5 to 0.6 mag (Peng et al. 2006). There- 
fore, there will be a color shift towards the red by about 
0.10-0.15 for a blue GC by the time that « 20-25% of 
the GC light originates from field stars. This effect will 
be notable mainly for the blue GC population, since the 
color of red GCs is much closer to that of field stars. 
As will become clear from the subsequent discussion, the 
capture rate per unit GC mass is higher for more mas- 
sive GCs, which will lead to a slope in color-magnitude 
space. This qualitatively explains finding (1) from above. 
In addition, the effect will be stronger in the inner part 
of the host galaxy because of the higher field star den- 
sity in such regions, consistent with finding (2). Thus, 
field star capture seems qualitatively consistent with the 
observed trend. 

We now quantify this scenario by considering a GC 
with a mass of 10 6 Mq, a mass that roughly corresponds 
to the brightest clusters in Groups 1 and 2 in Figure |2J 
For an assumed M/L=2.5, this corresponds to a total of 
4xl0 5 solar luminosities. We assume the GC to be on a 
circular orbit around M87. This is certain ly an extreme 
assumption, given that Cote et al. f2Q0\lJ) rule out such 
tangentially biased orbits for the GC system as a whole. 
We nevertheless use it for the sake of this simple estimate. 
We adopt a radius of 3 kpc, representative for the sample 
within 5 kpc projected galactocentric distance. 

We may estimate the number of captured stars by mul- 
tiplying the total number of stars within the volume of 
influence of the GC with the fraction of stars having rel- 
ative velocities below the GC escape velocity. As the 
volume of influence, we define a torus around the orbit 
of the GC with an inner radius of the GCs tidal ra- 
dius with respect to its host galaxy M87. For calculating 
the tidal radius, we use the ma ss pro file of M87 derived 
by Romanowsky & Kochanek (2001). They quote the 
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integrated mass at a distance of 3 kpc to be 1.1 xlO 11 
Mq, w hich i s comparable to the value found by Cote 
et al. i|200H) . Together with the 10 6 M Q mass of the 
GC, this gives a tidal radius of 9 pc. The total volume 
of the torus then is 4.8xl0 6 pc 3 . We use the same paper 
by Romanowsky & Kochanek to get the 3-dimensional 
stellar density of M87 at 3 kpc, which is about 0.06 L Q 
pc -3 . That is, along one orbit there will be a total of 2.9 
x 10 5 solar luminosities within the volume of influence of 
our test GC. 

However, only a small fraction of those will be cap- 
tured: the average escape velocity v esc of the GC within 
its tidal radius of 9pc is about 37 km s _1 , while the 
stellar velocity dispersion in M87 is about 275 km s _1 
(Romanowsky & Kochanek 2001). The GC population 
in M87 has a velocity dispersion of ~ 400 km s^ 1 . The 
probability distribution for the relative velocity between 
the field stars and a GC is given by the convolution of 
both distributions, which is a Gaussian of width a = 485 
km s _1 . The fraction of stars (in solar luminosities) that 
can be captured in the ID case is given by the area below 
the Gaussian between +/- which is 0.035. For the 

full 3D case, this is 0.035 3 = 4.3 x 10~ 5 , or one out of 
2.3 xlO 4 stars. Given that in total 2.9 xlO 5 stars are in 
the volume of influence in each orbit, this implies a cap- 
ture rate of one dozen stars per orbit 11 . For a Keplerian 
orbit at 3 kpc, one revolution lasts about 2.5 xlO 7 years. 
In a Hubble time, one can therefore capture about 6000 
stars. This is 1.5% of the total assumed GC luminosity, 
not sufficient to fully explain the observed effects. Only 
when considering a more extreme orbital distance of 1 
kpc does the total number of captured stars increase to 
about 4xl0 4 , i.e. ~10% of all stars. 

We now briefly examine the opposite effect of field star 
capture: namely cluster star evaporation. As is well 
known, mass segregation in GCs (Gunn & Griffin 1979; 
Jordan 2004) will result in evaporation affecting mainly 
the lowest mass stars. Lamers et al. (2006;) find for clus- 
ters of solar metallicity that evaporation will cause the 
overall GC color to become bluer with time (over most 
of their lifetime) due to the preferential loss of red, low- 
mass stars. Could such a color change help to explain 
our observations? Two effects are important: (1) the 
maximum blueward shift due to eva poratio n is larger for 
more massive clusters (Lamers et al. l)2006|) ). reaching up 
to A(V — I) ~ 0.10 mag; and (2) the dissolution time, 
tdiss, is longer for more massive clusters. To qualitatively 
explain our observations one would require the massive 
GCs to be in an early stage of evolution, so that evapora- 
tion has not yet had a significant effect. The fainter GCs 
would mainly define the color magnitude trend. How- 
ever, Figure |5] shows that the slope changes between in- 
ner and outer GCs are mainly driven by the brightest 
GCs. That is, evaporation is unlikely to account fully for 
the observed color-magnitude trend. Even more serious, 
evaporation cannot naturally explain the simultaneous 
presence of a trend for the blue GCs and its absence for 
the red. 

11 Note that close encounters with GC stars may indeed alter 
the capture rate, but an accurate assessment of this effect would 
require dedicated TV-body simulations, which is beyond the scope 
of this paper. 



5.5. GC Self Enrichment 

Generally speaking, self enrichment will occur in any 
stellar system, GCs included, if the system is able to 
retain a fraction of the gas expelled by the first genera- 
tion of supernovae and subsequently form a new genera- 
tion of stars. To create a metallicity-mass (or, equiv- 
alently, a color-magnitude) trend via self-enrichment, 
one would therefore require that more massive GCs are 
more efficiently self-enriched. The extent to which GCs 
in the Milky Way may have self-enriched, how ever, i s 
a contr oversi al subject: e.g., F rank fe Gisler (.1976), 
Smith 119961). Gnedin et al. il2002D. P armentier & 
Gilm ore (120011) . Dopita fc Sm ith jl98(| . Morgan & 
Lake l|1989j) . Thoul et al. l|200% 

In this context, it is interesting to note that the most 
massive GC in the Galaxy, uj Cen, is known to have 
a metallicity distribution function that shows a long 
tail towards high metallicities (e.g., Norris, Freeman & 
Mighell 1996; Stanford et al. 2006). Various stud- 
ies have found this high-metallicity population to be 
either coeval, or perhaps 2-4 Gyr younger, than the 
cluster's primary (metal-poor) population (e.g., Ferraro 
et al. 2004; Sollima et al. 2005; Hughes & Wallen- 
stein 2000; Hilker & Richtler 2000; Rey et al. 2004, Stan- 
ford et al. 2006). This is consistent with u> Cen having 
been at least partially self-enriched. Note, however, that 
because of its peculiar star formation history, numerous 
investigators have suggested that the cluster may be the 
stripped nucleus of a now-dissolved dwarf galaxy (e.g., 
Freeman 1993; Hilker & Richtler 2000; Bekki & Free- 
man 2003). If true, then it may be incorrect to consider 
ui Cen to be a typical example of a self enriched GC. 

This caveat aside, if self enrichment is common to the 
formation of the most massive GCs, then it may indeed 
help to explain our observations. Most significantly, the 
lack of a correlation among the red GC subpopulation 
is naturally explained in this scenario: i.e., the redder 
GCs are already metal-enriched, and a fixed amount 
dZ of supernovae ejecta will cause a relatively weaker 
metallicity increase, c?[Fe/H], for metal-rich GCs. For 
instance, for a system with [Fe/H] = -2 (Z = 0.0002), 
adding a fixed amount of metals dZ = 0.0002 will result 
in d[Fe/H]=+0.3 dex. According to Peng et al. (2006) 
this corresponds to a shift in (g — z) of about 0.06 mag 
and is about the observed amount of the color magnitude 
trend. However, for a system with [Fc/H] = —0.5 and 
the same dZ, we only have G?[Fe/H]=+0.013 dex, with 
a corresponding color change in {g — z) well below 0.01 
mag. 

A possible shortcoming of this self-enrichment scenario 
is that is does not explicitly account for environmen- 
tal effects such as a dependence on galactocentric ra- 
dius. On the other hand, feedback induced by pressure- 
confinement from the surrounding medium may be able 
to help reconcile this scenario with observations. 

5.6. The Distribution of GCs in Metallicity- Mass Space 

To conclude this discussion, we show in Figure El a 
metallicity-luminosity plot for the co-added sample of 
GCs in Group 1. To estimate GC metallicities, we have 
applied a quartic color-metallicity transformatio n base d 
on the empirical calibration data from Peng et al. (2006) . 
Note that the resulting metallicity distribution is not bi- 
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Fig. 12. — Metallicity plotted as a function of stellar mass 
and M z for GCs in Group 1, assuming a quartic color-metallicity 
transformation based on the empirical calibration data of Peng 
et al. 12006). Conv ersi on from luminosity to stellar mass is per- 
formed as in Figure ITol The vertical dashed line corresponds to 
the limiting color between blue and red GC subpopulations from 
Figure [S] The horizontal dashed line corresponds to the dividing 
point between brighter and fainter GCs from Figure [5] A KS test 
shows that the metallicity distributions in quadrant A and C stem 
from the same underlying distribution with only 0.1% probability. 
This probability is 91% for the color distributions in quadrants B 
and D. The probability for the joint sample of quadrants A and B 
compared to the joint samples C and D is 16%. 

modal, but rather shows a single pronounced peak at 
[Fe/H]~ —0.4 dex and a long, extended tail towards 
lower metallicities. This confirms recent warnings (e.g., 
Richtler 2005; Yoon, Yi & Lee 2006) that bimodality in 
color does not necessarily imply bimodality in metallic- 
ity. 

Using the data in Figure El we can compare the 
metallicity distributions of GCs brighter and fainter 
than M z = —10 mag. For the metal-poor popula- 
tion, ([Fe/H]<-0.6), the fainter GCs have a significantly 
broader tail towards lower metallicities than do the 
bright ones. A KS-test gives a 0.1% probability for a 
common underlying distribution between both samples. 
In the metal-rich part of the diagram, the agreement is 
very good, with a KS probability of 91%. In other words, 
the empirical finding of a color-magnitude trend for the 
blue GCs can hence be re-phrased in terms of metal- 
licity: the GC systems appear to have a skewed, but 
unimodal, metallicity distribution function whose metal- 
poor tail becomes increasingly extended at lower lumi- 
nosities/masses. 

6. SUMMARY AND CONCLUSIONS 

In this paper we have analysed the color-magnitude 
relation for the GC subpopulations in early-type galax- 
ies belonging to the Virgo Cluster. Our sample con- 
sists of 79 galaxies, spanning the range —21.7 < Mb < 
— 15.2 mag, drawn from the ACS Virgo Cluster Survey 
(Cote et al. 2004). Our principal findings can be sum- 
marized as follows: 

1. Based on KMM fits, we find a highly significant cor- 
relation 7z = = -0.037 ± 0.004 for the blue 



GCs in the co-added CMD of the three bright- 
est Virgo Cluster galaxies (M49, M87, M60). The 
sense of the correlation is such that brightest GCs 
in this subpopulation are redder than their faint 
counterparts. For M87 and M60 alone, we find 
comparable correlations; by contrast, the blue GCs 
in M49 do not appear to follow such a trend. In no 
galaxy do we find a significant correlation between 
(g — z) and M z for the GCs associated with the red 
subpopulation. 

2. Based on KMM fits, we find the slope j g = d[9 ~ g z) 
within the blue subpopulation to be much weaker 
than "f z . In turn, we find a mild positive correlation 
7 g for the red subpopulation, in the opposite sense 
to what is found for the blue subpopulation. 

3. We test KMM on artificial CMDs with, and with- 
out, implemented slopes for the blue subpopula- 
tion, and no slope for the red population. KMM cor- 
rectly reproduces the implemented slopes when us- 
ing the CMD (g - z) vs. M z . For (g - z) vs. M g , 
KMM gives a biased estimate for the color of the blue 
population and yields a positive slope for the red 
subpopulation, although no slope had been imple- 
mented. We trace this bias to the fact that in M g - 
unlike in M z - the luminosity function of blue GCs 
is populated to significantly brighter luminosities 
than for the red GCs. We conclude that for the 
real data, the strong slope j z is not an artifact of 
the fitting method, while the weaker slope j g is an 
artificial feature of the KMM fit. 

4. We derive correlations independent of KMM by us- 
ing a biweight estimator for GCs divided on the 
basis of a luminosity independent color. The re- 
sults confirm the significant slope found by KMM in 
M z and reveal no significant discrepancies between 
the slope measured using M z and M g . 

5. The slope measured from the blue subpopulation 
decreases with host galaxy luminosity. However, 
even among the faintest galaxies (—18.4 < Mg < 
— 15.2 mag) the slope is nonzero. 

6. The color-magnitude relation exhibited by the blue 
GCs is stronger for clusters with smaller projected 
galactocentric distances. This difference is driven 
mainly by GCs with M z < — 10 mag. 

We examine a number of physical mechanisms that 
might give rise to a color-magnitude relation with these 
characteristics: (1) contamination by super-star-clusters, 
stripped galactic nuclei, or ultra-compact dwarf galax- 
ies; (2) accretion of GCs from low-mass galaxies; (3) 
stochastic effects; (4) capture of galactic field stars by 
GCs; and (5) self-enrichment of individual GCs. Neither 
contaminants nor accretion of low-mass galaxies appear 
likely to fully explain our observations. Self-enrichment 
in more massive GCs may be able to account for the ob- 
served color-magnitude relation, although more work is 
needed to determine if the dependence on galactocentric 
radius can be explained within the context of the self- 
enrichment scenario. Although the capture of field stars 
does lead to a dependence on galactocentric radius that 
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is at least qualitatively in agreement with our observa- 
tions, this mechanism does not appear efficient enough 
to explain the overall size of the observed trend. 

While none of these scenarios can fully account for 
the full set of observations in a straightforward manner, 
self-enrichment and capture of field stars seem the most 
promising approaches. There is, of course, no a priori 
reason that a single mechanism is responsible for the ob- 
served color-magnitude relation, and conceptually, the 
combined effects of self-enrichment and field star cap- 
ture may be able to explain the observations: i.e., field 
star capture could act as a second order effect, produc- 
ing environmental variations of a global mass-metallicity 
relation induced by self-enrichment. This would allevi- 
ate the discrepancy between real and required field star 
capture rate estimated in this work. Dedicated N-body 
simulations going beyond these simple estimates will be 
required to accurately address the efficiency of field star 
capture in different environments. 
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